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Abstract 

The binding energy of three identical bosons is estimated by coupled differential 
equations which generalise the Feshbach-Rubinow approximation. This method 
turns out to be rather efficient, especially in the limit of vanishing binding. 
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1. Introduction 



There is no shortage of methods for solving the quantum- mechanical three-body 
problem: hyperspherical expansion, Faddeev equations, variational methods with 
sophisticated search for the free parameters, Monte-Carlo algorithms, etc. Each 
of them has merits and limitations. For instance, the hyperspherical expansion 
is a systematic scan of the Hilbert space, and can come in several variants such 
as "potential harmonics" or "adiabatic". However, it involves a rather delicate 
evaluation of the matrix elements of the potential within the basis of harmonics. 

In the 50's, a cleaver approximation, hereafter referred to as FR 1? was pro- 
posed by Feshbach and Rubinow [1]. It consists of seeking the bound-state wave 
function of three identical bosons of mass m as 

u(r) x + y + z 
*( W ) = ^, r= |— , (1) 

where x = |r 2 — r 3 |, y = |r 3 — ri| and z = |ri — r 2 | measure the distances between 
particles. The S-wave ground state, indeed, depends on three scalar, translation- 
invariant, variables, which can be chosen as x, y and z. The approximation consists 
of restricting to a function of their sum. Using astutely the stationnary properties 
of the Schrodinger equation, Feshbach and Rubinow derived the radial function 
obeyed by the radial function u(r) and energy E, 

15 14m 

-u"{r) + p«(r) + 2mV e (r)u(r) = — u {r)E , (2) 

where the effective potential results from the projection of the interaction V, 

V e (r) = I V(x,y,z)dr , j dr = -g / ydy (2r-y-z)zdz. (3) 

J J T JQ Jr—y 



This method has been applied to a variety of problems, in particular in nu- 
clear physics [2], and is still used occasionally. However, in the case of baryon 
spectroscopy [3], the energy provided by FRi was found [4] not to be very accu- 
rate, as compared, e.g., to the hyperspherical expansion in the approximation of 
a single partial wave, the so-called hyperscalar approximation (HA) . 

The best convincing success of FRi is perhaps in the domain of loosely-bound 
systems as encountered for instance in studying Borromean binding [5] . The FRi 
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approximation gives a good estimate of the coupling threshold" (73, to be compared 
to the coupling threshold gi of two-body systems. 

In Fig. 1, a comparison is made of the three-body binding energy of the simple 
pairwise exponential potential V{x,y,z) = gJ2 x v ( x )' with v(x) = — exp(— x), 
computed from FRi and HA, as a function of the coupling g. The constituent 
mass is set to m = 1. The FR4 method is clearly much better when the coupling 
threshold is approached, at g ~ 1.2. 
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Fig. 1: Ratio of Feshbach-Rubinow to hyperscalar ap- 
proximations of the three-body energy in a pairwise po- 
tential —g J2i<j ex P( — r ij)j as a function of the coupling g. 

The FRi method has been improved by Nogami et al. [2], Rosenthal et al. 
[6], and in particular adapted to deal with non-identical particles. In spite of its 
simplicity, the Feshbach-Rubinow method was never extensively used, probably 
because unlike the HA and other standard methods, it was never presented as 
the starting point of a converging expansion. We shall suggest below a possible 
remedy where the three-body wave function is written on a basis of orthogonal 
polynomials Pi with weight functions u n (r) that generalise the single u{r) of FRi 
and obey coupled equations. 

a 3>33 ensures a three-body bound state in the potential y] . v(\ri—rj\, where v is 
attractive or at least, contains attractive parts 

3 




The ground-state energy obtained from these coupled equations will be com- 
pared to the value given by standard variational methods. For estimating the 
quality of the wave function, the short-range correlation coefficient (^^(ri — r2)) 
will also be calculated. This quantity enters a number of decay or production 
rates, and is notoriously hard to compute accurately. Is estimate is sometimes 
made easier by using identities involving the derivative of the potential and cen- 
trifugal terms [4]. We shall restrict ourselves in this paper to a direct reading of 
the wave function at ri = r 2 . 



2. Formalism 

The FR4 ansatz is generalised as FR n 

nx,y,z) = ^P 1 (x,y,z) + ^P 2 (x,y,z) + --- , (4) 

where P n is a real and homogeneous polynomial of degree n which is symmetrical. 
The first one is simply Pi = r, and a truncation there corresponds to the original 
FRi method. Among the polynomials of degree n = 1 in x, y and z, only Pi oc 
x + y + z, indeed, is symmetric. For n > 1, the polynomials r n_1 Pi, . . . , rP n _i are 
already included in the first terms of the expansion (4), hence only a few genuinely 
new polynomials have to be introduced. b It is convenient to arrange the P n to be 
orthogonal. For n = 2, the new polynomial is 

P 2 ( x , y ,z) = ^-[72r 2 -49(x 2 + y 2 + z 2 )] , (5) 

where the coefficient is chosen such that J drP 2 = r 2 J drP 2 . 

If the expansion is truncated there, for testing purposes, the following set of 
coupled equations can be derived by straightforward calculus: 



h H 2~ui + a 12 u 2 H 2~w 2 +712 h Vuut + V 12 u 2 = — < 

T T T 10 



7 ( 6 ) 

a 22 u 2 H T u 2 + a 21 u x H T u\ +721 h V 2 \U\ + V 22 u 2 = —eu 2 , 

r r z r 15 



^ For three variables x, y and 2, the number p of symmetric polynomials of degree n 
does not exceed n till n=5, and (p—n)/n remains small for n>5 
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with e = 2m£, g n = -1 , a 12 = a 2 i = -2/(5 v / 237), a 22 = -389/395, /3n = 15/4, 
P12 = -(93/10) v^/79, (3 21 = -95/(2 v / 237), 711 = 722 = 0, 7l2 = - 721 = 
—98/ (5a/237). The 2x2 effective potential V^- results from the projections 



The coupled equations can be solved by several methods, for instance the 
discretisation procedure described in [4]. In short, r e [0, +00 [ is mapped into 
x G [0, 1[, and if iti(r) = each is expanded in Fourier series Vi = 

J2 n C n sm(nirx). The coefficients C n are translated into the values Vi(x n ) of the 
radial fonctions at equally-spaced points x n = nj (N + 1), where iV is the number 
of points. The system (6) is reduced to a 2N x 2N matrix equation for the quan- 
tities Vi(x n ), and its eigenvalues give a good approximation to the energy levels of 




(7) 



(6). 



The coefficient of short-range correlation, 



5= (#|5 (3) (x|tf>/(^> , 



(8) 



is given by 




for FRi and FR 2 , respectively. 
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3. Results 

The method has been first tested on the harmonic oscillator, v(x) = x 2 for which 
the exact solution is E = 6^/3/2 ~ 7.3485. The FRi gives E 1 = (36/7)^15/7 ~ 
7.5284, and the FR 2 leads to a remarkable improvement, with E 2 = 7.3537. Mean- 
while, the ratio of short-range correlation 5 to the exact value 5 ex = (3/2) 3 / 4 7r -3 / 2 
evolves from Si/5 ex = 1.57 for FRi to S 2 /5 ex = 0.968, again a dramatic betterment. 

The linear potential Yl,i<i r ij /2 = r is diagonal in the basis of polynomials 
Pi, but is not exactly solvable because the kinetic energy is not diagonal. The 
results are E\ ~ 3.906 for FRi and E 2 — 3.8635 for FR 2 , to be compared with 
El = q ~ 3.8647 and El<a — 3.8633 in the hyperspherical expansion with one and 
two partial waves, respectively [4]. The correlation coefficient is estimated to be 
<5i = 0.08488 for FRi and 5 2 = 0.05592 for FR 2 . For comparison, 5 ~ 0.05689 is 
obtained from four coupled hyperspherical waves [4], and 5 ~ 0.05702 from 176 
correlated Gaussians, a method to be presented shortly. 

As for the short-range potentials, let us take again the example of the ex- 
ponential potential v(x) = —gexp(-x) already considered in the introduction. 
Figure 2 displays the ratio of FR 2 to FRi ground-state energy as a function of the 
coupling g: the improvement saturates at about 2% at large coupling and becomes 
more interesting in the domain of very weak binding. 
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Fig. 2: Ratio ofFR 2 to FR 1 ap- 
proximations of the three-body 
energy in a pairwise potential 
-9j2i<j ex P(- r ij), as a func- 
tion of the coupling g. 



In Fig. 3, our results are compared with these of one the fashionable methods 
on the market. It consists of expanding the wave function on a basis of correlated 
Gaussians 

^(x, ? /,^) = ^C l {exp[-(a i v 2 + 26 l v.w + c l w 2 )/2] +•••} , (10) 

i 
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written in terms of the Jacobi variables 



v = r 2 - ri , 



2r, 



w = 



- £1 - r 2 



(ii) 



where the dots are meant for Gaussians deduced by permutation of the particles. 
Varga and Suzuki [7] have pushed the method of stochastic search of the parame- 
ters to a high degree of efficiency, and obtained very accurate estimates of binding 
energies in a variety of domains. We shall use here the strategy of Hiyama et al. [8] : 
Gaussians with diagonal quadratic forms, i.e., a^v 2 + qw 2 , are first introduced, 
and if aj 7^ q, they are supplemented by the terms deduced by permutation, which 
includes an angular dependence through v.w and carry the same weight factor Cj 
to ensure Bose symmetry. The range parameters ai and q are chosen to be any 
member of a geometric series {a, ar, ar 2 , . . . , ar^ -1 }. This means that, however 
large is the number of terms, only two non-linear parameters, a and r, are to be 
optimised. Note that N = 2 range parameters correspond to N' = 8 Gaussians, 
N = 3 to N' = 21, N = 4 to N' = 40, i.e., a fairly large number of terms. 

The plots in Fig. 3 indicate that the FR2 value is better than this variational 
estimate with N = 2, especially at vanishing binding, comparable to N = 3, and 
is defeated shortly by N = 4. This means that FR 2 is roughly equivalent to an 
expansion involving 40 correlated Gaussians, at least for this particular potential. 
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Fig. 3: Comparison of the energy in the exponential potential, as a 
function of the coupling g, for the FR 2 approximation (solid line), and 
the Gaussian expansion (dashed line) with N = 2 range parameters 
(left), N = 3 (centre) and N = 4 (right), corresponding to N' = 
8, 20, 40 terms. For clarity, the energies are shown by their ratio to 
the FRi approximation. 
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For further checking, a basis of correlated exponential functions has also been 
used, namely 



*f>(x,y,z) = ^2Ci [exp( 



aiX - biy - Ciz) H ] , 



(12) 



where the dots mean terms deduced by circular permutations of {a^, bi, q}. Wave 
functions of this type have been used in pioneering works on few-charge systems 
and in recent high-precision calculations of atomic systems [9]. The strategy of 
Ref. [8] can also be applied, to avoid numerical instabilities in the minimisation 
process, by restricting all range parameters a^, bi, . . . to belong to a single geometric 
series. Note that with suitable constraints on these range parameters, this method 
can mimic the FR1 method (o^ = bi = Ci Vj) or its generalisation [6] by Rosenthal 
and Haracz (RH), ({aj, bi, q} oc {a\, b\, c{\ Vz > 1), r 5 / 2 -u(r) being described as a 
sum of exponentials. This provides a check of the numerical solution of the radial 
FR1 equation or its RH analogue. 



In Ref. [6], indeed, Rosenthal et al. proposed to replace the symmetric ansatz 
F(r) = u(r) /r 5 / 2 , 2r = x + y + z of FR1 by 



where the dots correspond to permutation of x, y and z to restore the symme- 
try of the wave function. This method results in an integro-differential equation 
for F, whose solution is further optimised by varying r)2 and r]s empirically (a 
normalisation rji = 1/2 can be adopted). 

In Table 1, a comparison is made of the values obtained for the ground-state 
energy and correlation coefficient in the exponential potential at coupling g = 1.4 
and of the critical coupling g 3 at which three-body binding occurs. FRi and FR 2 
are compared to RH, as well as to the Gaussian expansion GN with up to iV 
different values of the range parameters, and it analogue EN with exponentials. 

At least with this potential, and for this domain of coupling close to the sta- 
bility threshold, the generalisation of the Feshbach-Rubinow approximation with 
two coupled differential equations, FR2, appears as more efficient than the single 
integro-differential equation RH. The ground-state energy, short-range correlation 
and coupling threshold approaches the exact value to about 1% accuracy for the 
simple exponential potential. 



^(x, y, z) = F(rnx + i] 2 y + i] 3 z) H , 



(13) 
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Table 1. Ground-state energy E3 and correlation coefficient 5 for an 
exponential potential with coupling g = 1.4, and critical coupling gs 
for three-body binding, with FRi, FR2 and RH methods and variational 
calculations using correlated Gaussian (G ) or exponential (E ) basis func- 
tions with varying number of allowed range parameters. 





G2 


GA 


G8 


E2 


EA 


E5 


FRi 


RH 


FR 2 


—E 3 
5 

93 


0.03055 
0.00840 
1.2194 


0.03575 
0.00972 
1.1613 


0.03586 
0.01000 
1.1563 


0.03466 
0.01251 
1.2158 


0.03583 
0.01005 
1.1587 


0.03586 
0.01002 
1.1560 


0.03468 
0.01212 
1.1751 


0.03475 
0.01224 
1.1649 


0.03542 
0.01091 
1.1644 



4. Outlook 

In this paper, the Feshbach-Rubinow method is generalised by expanding the 
ground-state wave-function of three identical bosons on orthogonal polynomials 
of the interparticle distances. The coefficients of these polynomials are functions 
of the sum of distances and obey coupled equations which can be derived by 
calculus and solved numerically. The results obtained with two equations provide a 
dramatic improvement in the case of confining potentials where the single-equation 
version of this method was known to be inaccurate. For the short-range potentials, 
the results are also significantly better, and become really excellent in the limit of 
loose binding, 

Further investigations would be necessary to check the convergence when the 
number of equations is increased, and for excited states. The accuracy of the 
method should also be probed for potentials with hard core, which often occur in 
applications. The case of particles with different masses or interaction properties 
also deserves further work. The binding energy and wave function of the positro- 
nium ion and other two-electron atoms, which were the subject of interesting 
developments of the Feshbach-Rubinow method, will be studied in a forthcoming 
paper within the formalism of coupled equations. 

Acknowledgments: I would like to thank Elian Masnada for stimulating discus- 
sions that led to this study, and Muhammad Asghar for very useful comments on 
the manuscript. 
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